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HIGHLIGHTS 


•  An  analytical  model  capable  of  explaining  experiments  is  derived  in  a  natural  way. 

•  Mechanism  and  limiting  factor  for  Ni  particle  growth  are  discussed  and  quantified. 

•  Only  one  adjustable  parameter  describing  electrode  morphology  is  used  in  the  model. 

•  The  model  is  in  very  good  agreements  with  the  available  experimental  data  sets. 
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The  agglomeration  of  Ni  particles  in  nickel— yttria  stabilized  zirconia  (YSZ)  anode  is  an  important 
degradation  mechanism  for  the  solid  oxide  fuel  cell  and  is  widely  believed  to  be  driven  by  surface 
diffusion.  This  work  aims  to  develop  a  quantitative  model  to  describe  the  agglomeration  kinetics.  The 
model  treats  the  anode  as  a  system  of  random  packing  Ni  and  YSZ  particles.  Surface  diffusion  occurs 
between  the  connected  Ni  particles  of  different  sizes  characterized  by  two  representative  radii,  but  is 
influenced  by  the  YSZ  network.  The  Fields  law  for  diffusion,  the  Gibbs-Thomson  relation  for  vacancy 
concentration  and  the  coordination  number  theory  for  percolating  Ni  network  are  employed  in  the 
mathematical  derivation.  The  growth  kinetics  is  expressed  as  an  analytical  function  consisting  of  two 
model  parameters,  one  for  the  Ni-particle  size  distribution  and  the  other  for  the  influence  of  the  YSZ 
backbone.  The  model  is  in  excellent  agreement  with  the  available  experiments.  The  influence  of  the  YSZ 
backbone  is  further  considered  to  obtain  a  model  with  just  one  fitting  parameter.  The  one-parameter 
model  is  also  in  good  agreement  with  the  experiments  and  the  fundamental  physics  for  the  Ni- 
particle  growth  is  therefore  believed  to  be  well  characterized. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Reliability  and  long-term  stability  are  crucial  for  solid  oxide  fuel 
cell  (SOFC)  to  become  a  commercially  mature  energy  conversion 
technology.  Research  on  the  mechanisms  responsible  for  the  per¬ 
formance  degradation  is  receiving  increased  attention  in  recent 
years.  Degradation  occurs  for  all  SOFC  components  and  for  a  variety 
of  reasons.  For  example,  crack  formation  and  loss  of  ionic  conduc¬ 
tivity  in  the  electrolyte  [1—6],  increase  of  the  contact  resistance 
between  electrodes  and  current  collectors  [7-15],  reactions  be¬ 
tween  electrolyte  and  electrode  materials  and  densification  of  the 
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electrode  structure  [16,17],  coarsening  of  electrode  particles  [18- 
28]  ,  delamination  between  electrode  and  electrolyte  [7,29-31] 
and  the  Cr  poisoning  of  the  cathode  [32-34],  etc.  Among  the 
degradation  phenomena  observed,  the  agglomeration  of  Ni- 
particles  in  anode  of  SOFC  has  received  particular  attention  as  the 
porous  Ni-YSZ  cermet  is  the  most  widely  used  anode  material  and 
the  agglomeration  process  that  results  in  decreasing  electro¬ 
chemical  performance  may  occur  in  a  relatively  short  period  of  time 
[18-21,25-28].  Due  to  the  complexity  of  the  agglomeration  phe¬ 
nomena  and  their  underlying  mechanisms,  the  progress  in  the 
development  of  theoretical  models  to  describe  the  kinetics  of  this 
process  is  rather  limited. 

The  agglomeration  of  Ni-particles  is  strongly  affected  by  the 
contents  of  YSZ  and  water  vapor.  The  presence  of  YSZ  limits  the 
space  for  the  growth  of  the  Ni  particles  [21,26].  The  water  vapor 
accelerates  the  agglomeration  process,  though  its  effect  on  the  final 


S.  Gao  et  al.  /  Journal  of  Power  Sources  255  (2014)  144-150 


145 


grain  size  of  Ni  particles  may  be  limited  [8,35].  The  mechanism  for 
the  influence  of  water  vapor  on  the  agglomeration  process  is 
complicated  and  is  conjectured  to  be  due  to  its  role  in  changing  the 
partial  pressure  of  Ni  gas  species  35]  and  the  surface  diffusion 
coefficient  of  Ni  atoms  [36],  etc.  To  avoid  the  ambiguity  caused  by 
the  yet  to  be  understood  role  of  water  vapor,  the  existing  models 
focus  on  describing  the  growth  of  Ni  particles  working  under  low 
water  vapor  pressure.  A  general  feature  of  these  models 
[19,26,27,35,37,38]  is  to  augment  the  traditional  theory  for  the 
agglomeration  of  particles  with  some  empirical  function  to  account 
for  the  observed  asymptotic  behavior  of  Ni  particle  growth.  For 
example,  Va(3en  et  al.  [26]  developed  a  theory  for  the  growth  of  Ni 
particles  with  the  surface  diffusion  as  the  dominant  mechanism, 
but  had  to  arbitrarily  introduce  a  limiting  function  to  mimic  the 
asymptotic  feature  of  the  Ni  particle  growth.  Even  though  the 
theory  may  fit  well  with  their  own  experimental  data,  the  deduced 
surface  diffusion  coefficient  of  Ni  is  2—3  orders  of  magnitude 
smaller  than  the  literature  data.  Similar  limiting  function  was  also 
used  empirically  by  Nakajo  et  al.  [38]  in  their  modeling  of  the 
coarsening  of  Ni  particles.  Naturally,  it  is  desirable  to  avoid 
enforcing  the  experimentally  observed  asymptotic  feature  in  the 
theoretical  model  by  arbitrarily  introducing  some  empirical  func¬ 
tion.  Recently,  Chen  et  al.  applied  a  phase-field  approach  to  simu¬ 
late  the  microstructural  evolution  of  Ni-YSZ  anode  in  terms  of  the 
triple  phase  boundary  (TPB)  density  [27  .  Even  though  such  simu¬ 
lation  eliminates  the  use  of  the  empirical  limiting  function,  it  lacks 
the  transparence  possessed  by  an  analytical  model.  Moreover,  it  is 
nontrivial  to  compare  the  simulation  results  with  the  experimental 
observations.  For  example,  the  presented  simulation  results  for  the 
chosen  modeling  parameters  indicate  much  faster  agglomeration 
processes  than  typically  observed  experimentally. 

This  work  focuses  on  developing  a  coherent  theoretical  model 
for  the  agglomeration  of  Ni-particles  in  SOFC  anode.  The  model  in 
this  paper  builds  on  the  foundation  of  significant  earlier  literature. 
The  model  uses  the  surface  diffusion  as  the  dominant  agglomera¬ 
tion  mechanism  and  incorporates  the  concept  of  coordination 
number  in  the  percolation  theory  and  the  statistical  nature  of  the 
random  packing  of  Ni-YSZ  particles  in  the  anode.  The  resulting 
theory  is  in  a  simple  analytical  form  with  two  adjustable  parame¬ 
ters  and  fits  the  experiments  well.  The  two  parameters  are  further 
reduced  to  one  based  on  reasonable  physical  consideration.  The 
model  with  only  one  adjustable  parameter  is  also  shown  to  fit  the 
experiments  very  well. 


Fig.  1.  Schematic  of  the  agglomeration  of  Ni-particle:  a)  the  process  of  coarsening,  b) 
geometric  model  of  the  two-particle  system. 


thermodynamically-driven  diffusion  of  vacancy  from  the  large 
particle  to  the  small  one  (Fig.  lb).  Or  equivalently,  the  Ni  atoms 
diffuse  from  the  small  particle  to  the  large  one.  As  a  result,  the  size 
of  the  large  particle  grows  while  the  small  particle  gets  smaller. 

The  rate  of  volumetric  change  of  the  large  particle  should  be 
equal  to  the  total  volumetric  flow  diffused  from  nearby  particles: 

w =  47rr'2^  m  ^  H =  njvAsfi  (1) 


2.  Theory 

The  agglomeration  of  Ni  particles  in  Ni-YSZ  cermet  is  in  prin¬ 
ciple  a  capillarity-driven  phenomenon.  The  growth  of  Ni  grains  may 
take  place  through  different  mechanisms:  evaporation— conden¬ 
sation,  bulk  diffusion  and  surface  diffusion.  The  evaporation- 
condensation  mechanism  is  negligible  as  the  vapor  pressure  of  Ni- 
containing  species  is  very  low  at  the  typical  SOFC  operating  con¬ 
ditions  (T  <  1000  °C)  [35].  The  grain  boundary  or  bulk  diffusion  is 
estimated  to  be  3  orders  of  magnitude  smaller  than  the  surface 
diffusion  [27].  Therefore,  it  is  reasonable  to  consider  the  surface 
diffusion  as  the  only  driving  mechanism. 

The  mass  transport  via  the  surface  diffusion  takes  place  between 
Ni  particles  of  different  curvatures.  The  larger  Ni  particles  grow  at 
the  expense  of  the  smaller  ones  which  may  tend  to  disappear,  as 
schematically  illustrated  in  Fig.  la.  To  capture  the  essence  of  the 
physical  process  in  simple  terms,  the  complicated  geometry  of  Ni 
particles  in  a  Ni-YSZ  cermet  is  modeled  as  a  system  of  two  sets  of 
spherical  Ni  particles,  with  radii  of  r\  and  rs  for  the  large  and  small 
particles,  respectively.  Due  to  the  difference  in  the  vacancy  con¬ 
centrations  of  the  large  and  small  particles,  there  is  a 


where  V\  is  the  volume  of  the  large  Ni  particle,  t  is  the  time,  Vf  is  the 
volumetric  flow  between  the  large  and  small  Ni  particles,  n  is  the 
number  of  Ni  particles  connected  to  the  large  particle  for  the 
diffusion  to  take  place,  jv  is  the  magnitude  of  vacancy  flow  between 
a  contacting  pair  of  large  and  small  particles,  AS  is  the  surface 
diffusion  area  of  the  vacancy  flow,  Q  is  the  diffusion  volume 
element  that  is  often  reasonably  approximated  by  the  atomic  vol- 
ume  of  Ni  (10.9  A3). 

For  the  convenience  of  evaluating  the  quantities  in  Eq.  (1)  and 
connecting  the  theoretical  model  with  the  experiment,  it  is  helpful 
to  use  the  mean  particle  radius,  rm  =  (q  +  rs)/2,  and  introduce  a 
proportionality  parameter  (3  =  (q  -  rs)/2rm  regarding  the  particle 
radius  difference.  Following  Vapen  et  al.  [26],  (3  is  assumed  to  be 
constant  during  the  agglomeration  process.  The  assumption  may  be 
empirically  justified  by  the  fact  that  the  average  particle  radius 
difference  is  generally  large  when  the  average  particle  size  is  large. 
However,  it  should  be  pointed  out  that  (3  by  itself  is  only  a 
phenomenological  parameter  like  q  or  rs.  They  are  used  to  simplify 
the  kinetic  model  and  should  not  be  interpreted  too  literally.  The 
final  test  of  the  model  rests  on  the  comparison  of  the  theoretical 
and  experimental  results  on  the  Ni  particle  growth,  rm  =  rm(t). 
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In  terms  of  rm  and  /3,  one  has  q  =  (1  +  /3)rm  and  rs  =  (1  -  /3)rm. 
Eq.  (1 )  may  be  rewritten  as, 


7  yVNi/fVNi 

zNi,Ni  —  z  — / 

'  m  V '  m 


Vysz\ 

rYSZ  / 


(7) 


4u(l+/?)3r^  =  n/vASO.  (2) 

Based  on  a  random  packing  of  Ni-YSZ  spheres,  the  average 
number  of  small  particles  connected  to  one  large  particle  may  be 
expressed  as: 


n  = 


(3) 


where  ZNi  Ni  is  the  coordination  number  between  Ni  particles,  i.e., 
the  average  number  of  contacts  between  Ni  particles  39].  ZNi  Ni  is 
divided  by  2  because  only  the  contacted  small  particles  contribute 
to  the  diffusion  flow.  /<Y sz  is  used  to  account  for  the  blocking  effect 
of  the  YSZ  backbone  on  the  Ni  particle  growth. 

During  the  process  of  Ni  particle  growth,  the  YSZ  particles  are 
known  to  remain  essentially  intact  due  to  the  very  high  melting 
temperature  of  the  YSZ  material  [19-21  ].  The  backbone  formed  by 
the  unchanging  YSZ  particles  provides  a  limited  space  for  the  Ni 
particle  growth,  as  shown  schematically  in  Fig.  2.  Moreover,  the  YSZ 
backbone  also  affects  the  connectivity  of  the  percolating  network  of 
Ni  particles.  The  total  number  of  Ni  particles  is  expected  to  grad¬ 
ually  decrease  in  the  process  of  agglomeration.  During  the 
agglomeration,  some  initially  percolated  Ni  particles  may  become 
disconnected  from  the  percolating  network,  while  the  others 
remain  connected,  as  illustrated  in  Fig.  2.  The  connected  Ni  particles 
remain  active  in  the  agglomeration  process,  while  the  disconnected 
Ni  particles  stop  growing  in  size.  That  is,  there  is  some  probability 
for  a  particle  to  stop  growing  at  a  given  time.  In  other  words,  there 
is  a  distribution  of  growth  lifetime  for  the  Ni  particles  [40].  /<Y sz  is 
introduced  in  Eq.  (3)  to  describe  this  fact  and  corresponds  to  the 
probability  of  Ni  particle  with  a  growth  lifetime,  T,  larger  than  a 
given  time,  t,  i.e.,  kYSZ  =  P(T  >  t). 

If  the  probability  per  unit  time  of  a  Ni  particle  stopped  growing 
at  time  t  isA  that  may  be  reasonably  assumed  to  be  constant  [40], 
the  probability  of  Ni  particles  stopped  growing  in  a  time  interval  (t, 
t  +  At)  would  be: 


Zysz,ysz 


2  5ysz/ /^^n[  +  %szA 
rYSZ  \rm  fysz/ 


(8) 


where  ijNi  (%sz)  is  the  volume  fraction  of  Ni  (YSZ)  based  on  the 
solid  phase  alone  (r;Ni  +  %sz  =  1).  Z  is  the  average  coordination 
number  and  rYSZ  is  the  mean  YSZ  particle  radius.  As  YSZ  is  the 
majority  component  of  the  anode  materials  in  most  SOFC  designs 
and  the  YSZ  particles  are  known  to  experience  very  little  change 
over  time  [19-21  ],  it  is  reasonable  to  assume  that  the  coordination 
number  of  YSZ  particles,  ZYSZ)Y sz,  is  a  constant.  Combining  Eqs.  (7) 
and  (8)  then  gives, 


%i  ^Ni  u 

ZNi,Ni  =  zYSZ,YSZ,o^  =  =  —  (9) 

rY sz  rmo  'T  rYszo  m 

HereZ0  is  the  initial  average  coordination  number,  rm0  (rY SZo)  is 
the  initial  mean  radius  of  Ni  (YSZ)  particles  and  the  constant  B  is 
given  by: 


B  = 


Z°M 

f"m0 


VNi 


_l_  Vysz  ‘ 
rYszo 


(10) 


Hence,  the  time  evolution  of  the  average  number  of  small  par¬ 
ticles  connected  to  a  large  particle  (Eq.  (3))  is  given  by: 


n  =  e  (11) 

m 

The  vacancy  flow  between  Ni  particles  of  different  sizes  may  be 
determined  by  the  Fields  law: 


,  _  dC  AC 

Jv  -  Uv,s^~Uv,sM 


(12) 


wheredC  is  the  difference  in  the  vacancy  concentration,  Ax  stands 
for  the  diffusion  distance  and  DV)S  is  the  vacancy  surface  diffusion 
coefficient.  According  to  the  Gibbs-Thomson  relation  [42],  the 
vacancy  concentration  may  be  expressed  as: 


P(t  <  T  <  t  +  At|T  >  t)  =  A- At,  (At— >0).  (4) 

Here  the  condition  “T  >  t”  is  required  as  the  particles  stopped 
growing  during  (t,  t  +  At)  necessarily  means  that  their  lifetimes  are 
longer  than  T.  By  nature  of  the  conditional  probability,  Eq.  (4)  may 
be  rewritten  as: 


A  =  lim  P(t  <  T  <  t  +  AtlT  >  t) /At 

At— >0 

=  lim  P(t  <  T  <  t  + At)/P(T  >  t) /At 
_ lim  P(r>t)-P(T>t+At) 

-  At"  0  P(T>t)M 

_  1  i  m  feysz(t)— kysz(f+At) 

fcvsz(t)At 

=  -1  dfcysz(t) 

/<YSZ(t)  dt  • 

Solving  Eq.  (5)  yields: 


(5) 


/<YSZ (t)  =  e  (6) 

where  the  initial  condition  /cysz(t  =  0)  =  1  has  been  used. 

The  coordination  numbers  between  Ni  particles  and  YSZ  parti¬ 
cles  in  the  SOFC  anode  are  respectively  39,41  ], 


Cr  =  Ce  exp 


2yQ\ 

Wr 


=  Ce  1 


2yQl 
kBT  r 


(13) 


where  Ce  is  the  vacancy  concentration  for  a  flat  surface,  /<b  is  the 
Boltzmann  constant,  T  is  the  temperature,  7  is  the  surface  energy 
and  r  is  the  radius  of  Ni  particle.  The  difference  in  the  vacancy 
concentration  is  therefore  given  by: 


AC  =  Cn  -  Cr,  = 


ekRrU  r, 


(14) 


As  seen  in  Fig.  lb,  the  diffusion  distance  may  be  written  as: 


Fig.  2.  The  limitation  on  the  growth  space  of  Ni  particles  imposed  by  the  YSZ 
backbone. 
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Fig.  3.  Comparison  of  the  experimental  data  and  the  theoretical  fittings  (Eq.  (20)):  a)  Cell  A  [19],  b)  Cell  B  [20],  c)  Cell  C  [21]. 


toe  =  \/(r,  +  rs)2  +  (n  -  rs)2. 

Hence,  the  vacancy  flow  is  given  by: 


jv  =  Dv,sCf^ 


A  rm 


'CbT  (rm~^rm)  (rm+Arm) 


2yDs _ 

k°T  (l-/52)(l+<f)°^ 


(15) 


(16) 


where  Ds  =  Dv,sCe@  is  the  atomic  surface  diffusion  coefficient  [43]. 

The  surface  diffusion  area  for  the  vacancy  flow  may  be  roughly 
estimated  with  a  thickness  3S  of  the  inter-atomic  spacing  of  Ni 
(2.5  A)  and  a  circumference  2Tcrm: 

AS  =  2nrm6s.  (17) 


With  Eqs.  (11),  (16)  and  (17),  Eq.  (2)  may  be  rewritten  as: 


drm  e~Xt 

dt  r4 


(18) 


where, 


C  =  D, 


yQ5s 


2kuT 


1  - 


O2)(l+02)“<l+« 


7  VNi 
-Z0 


Vnl _ | _  Vysz  ' 

fmO  fYSZO 


(19) 


Solving  Eq.  (18)  gives  the  time  evolution  of  the  mean  particle 
radius  of  Ni  as: 


rm  -  {}  ~e  +  rm0 ^  •  (20) 

Eq.  (20)  is  satisfactorily  simple  as  it  involves  only  two  adjustable 
parameters,  X  and  C,  or  effectively  A  and  (3  as  other  quantities  in  C 
are  determined  by  the  anode  composition  and  the  properties  of  the 
Ni  particles. 


Before  moving  further  on  applying  the  model,  it  is  helpful  to  be 
reminded  with  a  brief  discussion  on  the  conditions  for  the  applica¬ 
bility  of  the  theory.  As  mentioned  above,  these  conditions  include:  (1 ) 
Fuel  with  low  water  vapor  content  to  avoid  the  complication  caused 
by  the  effect  of  water  vapor.  Several  experiments  have  shown  that  the 
water  vapor  may  significantly  speed  up  the  agglomeration  process  via 
increasing  the  surface  diffusion  coefficient  of  Ni,  but  does  not  change 
the  final  state  [8,35].  These  experimental  results  imply  the  possible 
applicability  of  the  model  for  the  cases  with  high  water  vapor  con¬ 
tents.  However,  that  will  require  assumptions  about  the  unknown 
effects  of  water  vapor  on  the  diffusion  process  and  is  not  pursued 
here.  Moreover,  it  is  prudent  not  to  assume  the  applicability  of  the 
model  to  the  cases  with  high  water  vapor  pressures  as  often 
encountered  in  practical  operating  cells.  Nevertheless,  the  model 
should  be  directly  applicable  to  cells  with  the  design  of  gradual  in¬ 
ternal  reforming  where  the  water  vapor  remains  at  a  low  level 
throughout  the  cells,  as  demonstrated  experimentally  [44]  and 
illustrated  theoretically  [45];  (2)  Anode  microstructure  describable  by 
the  random  packing  of  Ni  and  YSZ  particles.  Fortunately,  the  random 
packing  model  is  known  to  describe  most  Ni— YSZ  composite  anodes 
well,  as  repeatedly  demonstrated  in  its  long  history  of  successful 
applications  [39,46-51].  However,  revision  may  be  required  for 
application  to  some  new  anode  designs  such  as  the  nano-composite 
anode  with  the  core-shell  structure  [52]  and  the  Ni  nanoparticle 
infiltrated  anode  [53],  as  illustrated  by  their  new  electrical  conduc¬ 
tivity  properties  [51  -53];  (3)  The  formation  of  a  stable  YSZ  backbone. 
The  YSZ  content  should  be  sufficiently  high,  say,  over  35vol%,  to  form  a 
highly  percolated  YSZ  network  in  a  conventional  Ni-YSZ  composite 
electrode  [39].  If  the  YSZ  content  is  low,  a  significant  number  of  YSZ 
particles  isolated  by  surrounding  Ni  particles  may  change  their  lo¬ 
cations  in  the  process  of  Ni  agglomeration.  In  such  cases,  both  the  role 
of  YSZ  on  limiting  the  Ni  agglomeration  is  diminished  and  the  coor¬ 
dination  number  theory  used  in  the  mathematical  derivation  may  be 
invalid.  Fortunately,  almost  all  Ni-YSZ  anode  designs  meet  the  model 
requirement.  As  a  side  note,  it  may  be  worthy  pointing  out  that  YSZ 
may  in  principle  be  replaced  by  other  stable  ceramic  materials,  but  the 
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Table  1 

Experimental  conditions  and  literature  data  used  for  validating  the  theoretical 
model  for  the  Ni-particle  growth  in  SOFC  anodes.  The  model  parameters,  f$  and  A, 
obtained  by  fitting  are  also  shown. 


Name 

Cell  A 

Cell  B 

Cell  C 

Temperature  T  (I<)  [19-21] 

1123.15 

1273.15 

1073.15 

Volume  fraction  of 

0.39 

0.4 

0.39 

Ni  »|Ni  [19-21] 

Volume  fraction  of 

0.61 

0.6 

0.61 

YSZ  77YSZ  [19-21] 

Initial  radius  of  Ni  particle 

0.29 

1.02 

0.32 

rm0  (pm)  [19-21] 

Initial  radii  of  YSZ  particles 

0.265 

1.0 

0.3 

rYSZ0  (Pm)  [19-21] 

Atomic  surface  diffusion 

1.04  x  10~10 

6.35  x  10  10 

4.80  x  10  11 

coefficient  Ds  (m2  h  ^  [26] 

Volume  element  Q  (m3)  [26] 

1.09  x  10"29 

1.09  x  1029 

1.09  x  1029 

Surface  energy  y  (Jm-2)  [26] 

1.9 

1.9 

1.9 

Thickness  5S  (m)  [26] 

2.5  x  lO"10 

2.5  x  10"10 

2.5  x  10"10 

Average  initial  coordination 

6.7 

6.7 

6.7 

number  Z0  [50] 

Fitting  parameter  /? 

0.  83 

0.  94 

0.  75 

Fitting  parameter  A  (h-1) 

4.70  x  10~3 

3.02  x  10~4 * * * * * 

1.17  x  10  3 

C  (m5  h-1)  (evaluated 

4.21  x  10"36 

2.03  x  1034 

1.73  x  10  36 

with  Eq.  (19)) 


replacement  by  volatile  materials  such  CGO  that  undergoes 
morphological  change  would  alter  the  agglomeration  kinetics  [54] 
and  would  invalidate  the  model. 

3.  Results  and  discussion 


3.2.  Reduction  of  model  parameters 


Eq.  (20)  may  be  rewritten  using  a  new  pair  of  model  parameters 
(/^»  tmmax)  as: 

rm  =  [(rrnmax  -  rmo)  (l  -  exP(  -  KP,  rmmax)t)  +  ^mo)]  5 

(21) 

where  A(/?,rmmax)  =  5C(/3)/(rgimax  -  r*^)  with  C(£)  given  by  Eq. 
(19).  Alternatively,  we  have, 


^5  _  y5 

1  mmax  1  m 


r5  _  r 5 

1  mmax  '  m0 


exp 


5C 


-r; 


mO 


(22) 


Here  rmmax  is  the  maximal  rm  for  the  agglomeration  of  Ni  par¬ 
ticles.  The  growth  of  Ni  particle  is  due  to  the  merger  of  connected 
neighboring  particles.  It  is  plausible  to  estimate  rmmax  by  allowing 
the  examined  particle  to  absorb  all  of  its  connected  Ni-particle 
neighbors.  That  is,  rmmax  may  be  estimated  through  the  following 
relation: 


47T  o 
~Yr\o + 


M/  (  ^Ni 

/mO  v  mO 


3  Mmax’ 

(23) 


where  rlmax  =  (1  +  jd)rmmax.  Quite  often,  the  initial  Ni  and  YSZ 
particles  are  of  similar  sizes.  Simplifying  Eq.  (23)  with  rm0  =  rYsz 
and  some  simple  algebra  yields, 


3.1.  Comparison  with  experimental  results 

Due  to  the  great  effort  required  for  the  long  term  testing,  there 
are  only  a  few  experimental  data  sets  available  for  the  model  to 
compare  to.  For  the  convenience  of  referencing,  the  experimental 
results  by  Faes  et  al.  [19],  Simwonis  et  al.  [20]  and  Tanasini  et  al.  [21  ] 
are  referred  below  as  Cell  A,  Cell  B  and  Cell  C,  respectively. 

Fig.  3  shows  the  comparison  of  the  theoretical  model  with  the 
experimental  data.  The  property  parameters  required  by  the  model 
are  chosen  to  be  as  realistic  as  possible  and  summarized  in  Table  1. 
For  example,  the  surface  diffusion  coefficients  of  Ni  correspond  to 
the  newest  data  by  Jaunet  [26].  These  coefficients  were  obtained 
without  the  water  vapor  effect  and  might  moderately  underesti¬ 
mate  that  for  the  cells  tested  with  a  small  amount  (3-4%)  of  H2O. 
The  underestimation  may  give  rise  to  somewhat  overestimated  /3 
values,  as  seen  in  Eq.  (19).  As  shown  in  Fig.  3,  the  current  model 
with  reasonable  property  parameters  fits  the  experimental  data 
very  well.  The  present  model  is  clearly  superior  over  the  model  by 
Va(3en  et  al.  that  has  to  rely  on  arbitrarily  low  surface  diffusion 
coefficients  (over  two  orders  of  magnitude  smaller  than  the  Jaunet 
values)  to  reproduce  the  measured  low  growth  rates  [26]. 


Fig.  4.  Dependence  of  rmmax/rm 0  on  the  model  parameter  0 


(24) 


The  dependence  of  rmmax/rm0  on  0  for  the  typical  case  of 
rjNi  =  0.4  and  Z0  =  6.7  [50]  is  shown  in  Fig.  4.  As  may  be  seen  from 
Fig.  4,  rmmax/rm0  is  essentially  a  constant  for  0  >  0.5. 

With  Eq.  (24),  there  is  only  one  adjustable  parameter,  /3,  in  the 
theoretical  model.  For  convenience,  the  two-parameter  model  and 
the  one-parameter  model  are  named  as  Model  A  and  Model  B, 
respectively.  Fig.  5  shows  the  comparison  of  Model  B  with  the 
experimental  data.  It  is  natural  to  expect  that  the  fitting  quality  of 
Model  B  with  only  one  adjustable  parameter  is  lower  than  that  for 
Model  A.  As  shown  in  Fig.  5,  however,  Model  B  also  fits  the  ex¬ 
periments  very  well.  The  worst  fit  is  seen  for  Cell  A.  Nevertheless, 
the  largest  fitting  error  is  only  15  nm,  or  a  relative  error  of  4%,  well 
below  the  experimental  uncertainty.  As  it  is  rather  difficult  to 
reproduce  several  different  experimental  results  with  just  one 
adjustable  parameter,  the  good  fitting  quality  of  Model  B  is  highly 
satisfactory.  Moreover,  it  is  noted  that  the  fittings  by  Model  A  and 
Model  B  also  produce  reasonably  similar  results  for  the  unknown 
parameters.  For  example,  the  values  of  (3  determined  by  Model  A/ 
Model  B  are  0.83/0.76,  0.94/0.94  and  0.75/0.68  for  Cell  A,  Cell  B  and 
Cell  C,  respectively.  The  corresponding  values  in  pm  for  rmmax  are 
0.366/0.383,  1.35/1.35  and  0.404/0.423,  respectively.  Therefore, 
Model  B  is  a  good  reflection  of  Model  A.  It  is  recommendable  to  use 
Model  B  instead  of  Model  A  as  reducing  the  number  of  adjustable 
parameters  is  often  strongly  favorable  for  a  theoretical  model. 


4.  Summary 

A  theoretical  model  for  surface  diffusion  driven  Ni-particle 
agglomeration  in  SOFC  anode  is  proposed.  The  diffusion  process 
is  determined  by  the  Fields  law  and  the  Gibbs— Thomson  relation  for 
a  system  of  random  packing  spherical  Ni  particles  with  two 
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Fig.  5.  Fitting  of  the  experimental  data  by  Model  B:  a)  Cell  A  [19],  b)  Cell  B  [20],  c)  Cell  C  [21]. 


representative  radii.  A  lifetime  distribution  for  the  Ni-particle 
growth  by  the  surface  diffusion  is  introduced  to  account  for  the 
statistical  connectivity  nature  of  the  percolating  Ni-particle 
network  confined  by  the  YSZ  backbone.  The  mass  diffusion  is 
further  quantified  through  the  concept  of  coordination  number  in 
the  percolation  theory.  A  simple  analytical  expression  with  two 
adjustable  parameters  is  obtained.  The  model  is  shown  to  be  in 
excellent  agreements  with  the  experimental  results  for  different 
microstructure  compositions  and  operating  temperatures. 
Furthermore,  the  role  of  the  YSZ  backbone  on  limiting  Ni-particle 
agglomeration  is  used  to  estimate  the  maximal  size  of  Ni  particle, 
yielding  a  model  with  only  one  adjustable  parameter.  The  one- 
parameter  model  is  also  found  to  agree  with  the  experiments 
very  well,  strongly  suggesting  that  the  fundamental  physics  for  the 
Ni-particle  growth  is  well  represented  in  the  theory. 
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